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Single crystals of EUCU 2 AS 2 and EuCu 2 Sb 2 were grown from CuAs and CuSb self-flux, respec¬ 
tively. The crystallographic, magnetic, thermal and electronic transport properties of the single 
crystals were investigated by room-temperature x-ray diffraction (XRD), magnetic susceptibility \ 
versus temperature T, isothermal magnetization M versus magnetic field H, specific heat C P (T) 
and electrical resistivity p(T) measurements. EUCU 2 AS 2 crystallizes in the body-centered tetrago¬ 
nal ThCr 2 Si 2 -type structure (space group /4/mmm), whereas EuCu 2 Sb 2 crystallizes in the related 
primitive tetragonal CaBe 2 Ge 2 -type structure (space group PA/nmm). The energy-dispersive x-ray 
spectroscopy and XRD data for the EuCu 2 Sb 2 crystals showed the presence of vacancies on the Cu 
sites, yielding the actual composition EuCui. 82 Sb 2 . The p(T) and C P (T) data reveal metallic char¬ 
acter for both EuCu 2 As 2 and EuCui.g 2 Sb 2 . Antiferromagnetic (AFM) ordering is indicated from 
the x(T)i C p (T), and p(T) data for both EUCU 2 AS 2 (Tn = 17.5 K) and EuCui.g 2 Sb 2 (Tn = 5.1 K). 

In EuCui. 8 2 Sb 2 , the ordered-state x(T) and M(H) data suggest either a collinear A-type AFM 
ordering of Eu +2 spins S = 7/2 or a planar noncollinear AFM structure, with the ordered moments 
oriented in the tetragonal ab plane in either case. This ordered-moment orientation for the A-type 
AFM is consistent with calculations with magnetic dipole interactions. The anisotropic y(T) and 
isothermal M(H ) data for EUCU 2 AS 2 , also containing Eu +2 spins S = 7/2, strongly deviate from 
the predictions of molecular field theory for collinear AFM ordering and the AFM structure appears 
to be both noncollinear and noncoplanar. 

PACS numbers: 74.70.Xa, 75.50.Ee, 65.40.Ba, 72.15.Eb 


I. INTRODUCTION 

The observation of high-T c superconductivity in FeAs- 
based 122-type compounds AFe 2 As 2 (A = Ca, Sr, Ba) af¬ 
ter suppressing the spin-density wave (SDW) transition is 
intriguing and these pnictides have been a topic of contin¬ 
uing research. 1-10 Because of their simple ThCr 2 Si 2 -type 
body-centered tetragonal crystal structure and the avail¬ 
ability of large single crystals, these compounds present 
a platform to understand the mechanism for the high-T/ 
superconductivity and other interesting properties of the 
iron-arsenide class of superconductors. The AFe 2 As 2 (A 
= Ca, Sr, Ba) compounds exhibit itinerant antiferromag¬ 
netic (AFM) spin-density wave (SDW) transitions that 
are accompanied by a tetragonal to orthorhombic struc¬ 
tural distortion, and superconductivity is realized by par¬ 
tial substitution at the A, Fe and/or As sites or by ap¬ 
plication of external pressure upon suppressing the SDW 
and structural transitions. 1-6,8 For example, K-doping at 
the Ba-site in BaFe 2 As 2 , which exhibits structural and 
SDW transitions near 140 K, 11,12 suppresses these transi¬ 
tions and results in superconductivity with T c up to 38 K 
for x ~ 0.4 in Bai_ a; K 2; Fe 2 As 2 . 1 Cobalt substitutions at 
the Fe site of BaFe 2 As 2 result in superconductivity with 
a maximum T c ~ 25 K at the optimum doping concen¬ 
tration x ~ 0.06 in Ba(Fei_ a; Co a ;) 2 As 2 . 5,13,14 Isovalent 
P-doping at the As site in BaFe 2 As 2 also leads to su¬ 
perconductivity in BaFe 2 (Asi_ 3: P a: ) 2 with T c ss 30 K for 
x « 0.3. 15,16 Superconductivity in BaFe 2 As 2 can also be 
induced by application of external pressure. 17 

If nonmagnetic Ba +2 in BaFe 2 As 2 is completely re¬ 


placed by Eu +2 which carries a local moment with spin 
S = 7/2, an interesting situation occurs in which both 
itinerant conduction carrier magnetic moments and lo¬ 
calized moments are present. In EuFe 2 As 2 , the Eu +2 
moments order antiferromagnetically below 19 K with 
an A-type AFM structure and the itinerant current car¬ 
riers undergo a SDW transition at 190 K with an associ¬ 
ated structural distortion. 18-21 In the A-type structure, 
the Eu ordered moments are ferromagnetically aligned 
in each afe-plane layer with the ordered moments aligned 
in the ab plane, but where the moments in adjacent lay¬ 
ers along the c axis are antiferromagnetically aligned (see 
Fig. 2(b) below for the proposed A-type AFM structure 
of the Eu ordered moments in EuCu 2 Sb 2 which is the 
same as in EuFe 2 As 2 ). The A-type AFM structure is 
therefore somewhat unusual, because the ferromagnetic 
(FM) alignment within an ab plane often arises from 
dominant intraplane FM interactions, and the AFM then 
arises from weaker AFM interplane interactions. The 
Curie-Weiss law for the magnetic susceptibility x is 

x(T) = ^22-, (1) 

where C is the Curie constant and 9 p is the Weiss tem¬ 
perature. In the present context we consider a Heisen¬ 
berg spin lattice with Hamiltonian PL = Y^<ij> Jij^i' Sj, 
where the sum is over distinct {S 2 , Sj} spin pairs 
and the exchange interactions Jij are positive for AFM 
interactions. For a spin lattice containing identical 
crystallographically-equivalent spins such as the Eu +2 
spins S = 7/2 in EuFe 2 As 2 (and in EuCu 2 As 2 and 



2 


EuCui. 82 Sb 2 , see the following), the Weiss temperature is 
given in general by molecular field theory (MFT) as 22-24 


S{S + 1) 
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where &b is Boltzmann’s constant and the sum is over all 
neighboring spins j of a given central spin i with exchange 
interactions Jij , respectively. Thus if the dominant inter¬ 
actions are FM (negative), then 9 p is positive (ferromag¬ 
neticlike) as in EuFe 2 As 2 , 21 even though EuFe 2 As 2 is an 
antiferromagnet. In this regard, we note that the simi¬ 
lar compound EuCu 2 P 2 does order ferromagnetically. 25 
The Eu spins-7/2 in EuCo 2 As 2 were claimed to exhibit 
A-type AFM ordering below 39 K. 26 

Like BaFe 2 As 2 , superconductivity is observed in 
EuFe 2 As 2 after suppression of the SDW and the as¬ 
sociated structural transition. The interaction of the 
Eu moments with superconductivity in pure and doped 
EuFe 2 As 2 has been extensively studied. 27-33 Partial sub¬ 
stitution of Eu by K leads to superconductivity in 
Eui_ a; K a ;Fe 2 As 2 with T c as high as 33 K for x = 
0.5. 27,28,33 Re-entrant superconductivity is observed in 
EuFe 2 As 2 on application of hydrostatic pressure. 30 Like 
substitution of Co for Fe in BaFe 2 As 2 , Co substitution for 
Fe in EuFe 2 As 2 also leads to superconductivity but the 
superconductivity in Eu(Fei_ a; Co a ;) 2 As 2 is re-entrant, re¬ 
vealing an important role of the Eu magnetic moment. 31 
However, Ni substitution for Fe in EuFe 2 As 2 does not 
induce superconductivity down to 2 K (Ref. 34) which 
contrasts with the observation of superconductivity in Ni- 
doped BaFe 2 As 2 . 35,36 Interestingly, Ni-doped EuFe 2 As 2 
is reported to exhibit FM ordering of the Eu moments 
below 20 K, 34 which is not too surprising in view of 
the dominant FM interactions in EuFe 2 As 2 as discussed 
above. It was observed that partial subsitution of Ni for 
Fe in Euo. 5 Ko.sFe 2 As 2 leads to a revival of the SDW tran¬ 
sition in Euo. 5 Ko. 5 (Fe 1 _ x Ni a ,) 2 As 2 as the holes in the sys¬ 
tem created by K-doping at the Eu site are compensated 
by the doped electrons introduced by partially replacing 
Ni for Fe. 37 

Recently we reported the crystallographic and physi¬ 
cal properties of SrCu 2 As 2 and SrCu 2 Sb 2 which are sp- 
band metals. 38 The Cu atoms are in a nonmagnetic 3d 10 
electronic configuration with a formal Cu +1 oxidation 
state. Here we report the physical properties of magnetic 
analogues of these compounds, namely EuCu 2 As 2 and 
EuCu 2 Sb 2 . EuCu 2 As 2 crystallizes in the body-centered 
tetragonal (bet) ThCr 2 Si 2 -type structure (space group 
IA/mmm ) and EuCu 2 Sb 2 in the primitive tetragonal 
CaBe 2 Ge 2 -type ( PA/nmm ) structure. 39 Diinner et al. 39 
reported the occurrence of vacancies on the Cu sites 
in EuCu 2 As 2 as deduced from their refinement of x- 
ray diffraction (XRD) data, yielding an actual compo¬ 
sition EuCui. 760 ( 5 )As 2 . No such vacancies were reported 
for EuCu 2 Sb 2 . The Eu atoms are crystallographically 
equivalent in both compounds. Sengupta et al. 40,41 re¬ 
ported magnetic and other physical properties of poly¬ 


crystalline EuCu 2 As 2 . They found from y(T) measure¬ 
ments that polycrystalline EuCu 2 As 2 orders antiferro- 
magnetically below Tn = 15 K, and from x(T) and 151 Eu 
Mossbauer measurements they found that the Eu ions 
are divalent (S = 7/2). 40 From the temperature of the 
peak of the zero-field-cooled x versus T at fixed H they 
inferred an approximately linear decrease in Tn with in¬ 
creasing H , where T^(H = 0) = 15.5 K and Tn —> 0 at 
H = 1.81 T. 41 From electrical resistivity p(T) measur- 
ments under pressure, they also inferred that Tn strongly 
increases with increasing pressure, reaching 49 K at a 
pressure of 10.7 GPa. 41 The crystallographic and mag¬ 
netic properties of the related compounds EuPd 2 As 2 
(Ref. 42) and EuPd 2 Sb 2 (Ref. 43) have also been re¬ 
ported, where AFM ordering of the Eu spins 7/2 is also 
found in each compound. The detailed physical proper¬ 
ties of EuCu 2 Sb 2 have not been reported before to our 
knowledge. 

Herein, we report the growth of EuCu 2 As 2 and 
EuCui.8 2 Sb 2 single crystals and their crystallographic, 
magnetic, thermal and electronic transport properties in¬ 
vestigated using powder XRD, x(T), magnetization M 
versus applied magnetic field H isotherms, specific heat 
C P (T), and electrical resistivity p(T) measurements. The 
experimental details are given in Sec. II and the crystal¬ 
lographic results in Sec. III. The physical property mea¬ 
surements of EuCui. 82 Sb 2 and EuCu 2 As 2 are presented 
and analyzed in Secs. IV and V, respectively. A summary 
is given in Sec. VI. 


II. EXPERIMENTAL DETAILS 

Single crystals of EuCu 2 As 2 and nominal EuCu 2 Sb 2 
were grown by the self-flux method. The high-purity el¬ 
ements Eu (Ames Lab), Cu 99.999% (Alfa Aesar), As 
99.99999% (Alfa Aesar) and Sb 99.999% (Alfa Aesar) 
were used for synthesis. Eu and prereacted flux (CuSb or 
CuAs) taken in a 1:5 molar ratio were placed in alumina 
crucibles which were then sealed in evacuated quartz 
tubes. The single crystals were grown by heating the 
quartz ampoules to 850 °C at a rate of 60 °C/h, held for 
5 h, heated to 1100 °C at a rate of 60 °C/h, held for 25 h 
and then slowly cooled to 800 °C at a rate of 2.5 °C/h. 
The crystals were separated from the flux by centrifuga¬ 
tion at that temperature. We obtained shiny plate-like 
crystals of typical size 2.5 x 2 x 0.3 mm 3 for EuCu 2 As 2 
and 4 x 3 x 0.5 mm 3 for EuCu 2 Sb 2 . 

The phase purity and the crystal structure of the 
crystals were determined by powder XRD collected on 
crushed single crystals using Cu Ka radiation on a 
Rigaku Geigerflex x-ray diffractometer. The single¬ 
phase nature of the crystals was further checked using 
a JEOL scanning electron microscope (SEM) equipped 
with an energy-dispersive x-ray spectroscopy (EDS) ana¬ 
lyzer. The single-phase nature of the crystals was inferred 
from the high resolution SEM images. Wavelength- 
dispersive x-ray spectroscopy (WDS) measurements were 



3 


also carried out to determine the compositions of the 
crystals. The average compositions obtained from the 
EDS and WDS measurements showed the expected 1:2:2 
stoichiometry for EUCU 2 AS 2 and EuCu 2 Sb 2 except for 
~ 9% vacancies on the Cu sites in EuCu 2 Sb 2 correspond¬ 
ing to an actual composition EuCui.g 3 ( 5 )Sb 2 . This lat¬ 
ter result is consistent with the site occupancies deter¬ 
mined from Rietveld refinements of the XRD patterns 
in Sec. Ill where a composition EuCui.g 2 (i) As 2 is found. 
This composition is contrary to the full occupancies re¬ 
ported in Ref. 39. The difference likely arises from differ¬ 
ences in sample preparation conditions. All molar mag¬ 
netic and heat capacity measurements reported here for 
EuCui.g 2 Sb 2 are normalized to a mole of EuCui.g 2 Sb 2 
formula units. 

The magnetization measurements were performed us¬ 
ing a Quantum Design, Inc., superconducting quantum 
interference device magnetic properties measurement sys¬ 
tem (MPMS). The output of the MPMS software for the 
magnetic moment of a sample is in Gaussian cgs elec¬ 
tromagnetic units (emu), which in terms of fundamental 
quantites is given by 1 emu = 1 G cm 3 . The magnetic 
field unit is the Oe, with 1 Oe = 1 G, and we use the 
Tesla (T) as a unit of convenience where 1 T = 10 4 Oe. 

Due to the large Eu spins-7/2 and the resulting demag¬ 
netizing fields that can be especially large at low temper¬ 
atures, the magnetic field H 0a applied in the ct th direc¬ 
tion has been corrected for the demagnetizing field H^a-, 
yielding the internal field H lnt in cgs units obtained from 

Hint a — H() a 47rA^daA/ct,, (3) 

where the demagnetizing factor Nd a is defined as for 
the SI system of units for which 0 < N^ a < 1 and 
x IVda = 1. Thus the M(H) isotherms are presented 
as M versus H[ n t isotherms and the reported intrinsic 
Xa = M a /Hi n t data are obtained from the observed 
^.obso. _ A-fobs Q /h 0o data using 

v obs 

Xa “ 1 - 47rAd a x° bs ' 

Since our crystal shapes can be approximated by rect¬ 
angular prisms, the magnetometric N,i a values were ob¬ 
tained using Eq. (1) in Ref. 44. Our reported magneti¬ 
zations and magnetic susceptiblities are normalized to a 
mole of formula units (f.u.), whereas M a in Eq. (3) and 
X° bs in the denominator of Eq. (4) are normalized to unit 
volume with respective units of G and dimensionless, re¬ 
spectively. In cgs units the volume-normalized quantities 
are obtained from the ones normalized per mole of f.u. 
by dividing by the molar volume Em [cm 3 /mol] as given 
in Table I below. 

The heat capacity was measured by a relaxation tech¬ 
nique using a Quantum Design, Inc., physical proper¬ 
ties measurement system (PPMS). The electrical resis¬ 
tivity measurements were performed using the standard 
four-probe ac technique using the ac-transport option of 
the PPMS with the current in the ab plane. Annealed 



FIG. 1: (Colour online) Powder x-ray diffraction patterns 
of (a) EuCu 2 As2 and (b) EuCu 2 Sb 2 recorded at room tem¬ 
perature. The solid line through the experimental points 
in (a) is the Rietveld refinement profile calculated for the 
ThCr 2 Si 2 -type body-centered tetragonal structure (space 
group /4/mmm), and in (b) the Rietveld refinement pro¬ 
file for the CaBe 2 Ge 2 -type primitive tetragonal structure 
(space group PA/nmm ). The short vertical bars mark the 
Bragg peak positions. The lowermost curves represent the 
differences between the experimental and calculated intensi¬ 
ties. The unindexed peaks marked with stars correspond to 
peaks from the flux. The Miller indices (hkQof the strongest 
peaks in (a) and (b) are indicated. Whereas all (hkf) com¬ 
binations are allowed for the primitive tetragonal structure 
of EuCu 2 Sb 2 in (b), only indices with h + k + l = even 
are allowed for the body-centered-tetragonal structure of 
EuCu 2 As 2 in (a). 


platinum wire (50 jim diameter) electrical leads were at¬ 
tached to the crystals using silver epoxy. 


III. CRYSTALLOGRAPHY 

The powder XRD data were collected on crushed sin¬ 
gle crystals of EuCu 2 Sb 2 and EUCU 2 AS 2 at room tern- 
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FIG. 2: (Color online) (a) ThCr 2 Si 2 -type body-centered 
tetragonal crystal structure (74 /mmm) of EUCU 2 AS 2 . 
(b) CaBe 2 Ge 2 -type primitive tetragonal crystal structure 
(P4/nmm) of EuCu 2 Sb 2 - The order of the Cu and Sb lay¬ 
ers in the lower half of the unit cell is reversed in (b) with 
respect to (a). In order to compare the structures, the ori¬ 
gin of the EuCu 2 Sb 2 unit cell is shifted by (1/4, 1/4, 1/4) 
from the atomic coordinates in Table II. The arrows in (b) 
show the ordered Eu moments in an A-type AFM structure, 
which is one of two AFM structures proposed in this paper 
for EuCui.8 2 Sb 2 . The other is a planar helix structure with 
the helix axis along the c axis. 


perature as shown in Fig. 1. Also shown are the struc¬ 
tural Rietveld refinement profiles of the XRD data using 
FullProf 45 software. The refinements of the XRD data 
indicate that the crystals are single-phase since no ex¬ 
tra peaks beyond those of the respective 122-type phases 
were observed except for the weak unindexed peaks 
marked with stars in the XRD pattern of EUCU 2 AS 2 
that arise from a small amount of adventitious flux on 
the surfaces of the crystals. Our Rietveld refinement 
confirmed the ThCr 2 Si 2 -type bet structure (74/mmm) 
for EUCU 2 AS 2 , and the CaBe 2 Ge 2 -type primitive tetrag¬ 
onal crystal structure (P4/nmm) for EuCu 2 Sb 2 - The 
two crystal structures are shown in Fig. 2, and are both 
ternary derivatives of the binary BaAU-type structure, 46 
consisting of layers of Eu, Cu and As/Sb atoms stacked 
along the tetragonal c axis. From Fig. 2, these two struc¬ 
tures differ in the arrangement of layers of Cu and As/Sb 
layers. The order of the Cu and Sb layers in the lower 
half of the unit cell is reversed in the CaBe 2 Ge 2 -type 
structure of EuCu 2 Sb 2 in Fig. 2(b) compared to that in 
ThCr 2 Si 2 -type EUCU 2 AS 2 in Fig. 2(a), thus resulting in 
a loss of the mirror plane perpendicular to the c axis and 
of the inversion symmetry about the unique Eu position 
in the CaBe 2 Ge 2 -type structure. While all Cu atoms oc¬ 
cupy equivalent Wyckoff 4 d positions in the ThCr 2 Si 2 
structure of EuCu 2 As 2 , in the CaBe 2 Ge 2 structure of 
EuCu 2 Sb 2 the Cu atoms are equally distributed between 
the 2 a and 2c positions. Thus in EuCu 2 Sb 2 there are two 
distinct types of Cu square lattices which have different 
lattice parameters and are rotated by 45° with respect to 


TABLE I: Crystallographic and Rietveld refinement param¬ 
eters obtained from powder XRD data for crushed crystals 
of body-centered tetragonal ThCr 2 Si 2 -type EuCu 2 As 2 (space 
group 74/mmm) and primitive tetragonal CaBe 2 Ge 2 -type 
EuCu 2 Sb 2 (space group PA/nmm). The molar volume Vm 
(volume per mole of formula units) is also listed. Shown for 
comparison are the lattice parameters for single crystals from 
Ref. 39 and for a poly crystalline sample of EUCU 2 AS 2 from 
Ref. 40. 



EuCu 2 As 2 



Lattice parameters 

This Work 

Ref. 39 

Ref. 40 

a (A) 

4.2330(1) 

4.215(1) 

4.260(1) 

c (A) 

10.1683(3) 

10.185(2) 

10.203(1) 

c/a 

2.4022(1) 

2.416(1) 

2.395(1) 

Well (A 3 ) 

182.20(1) 

180.95(12) 

185.16(11) 

Vm (cm 3 /mol) 

54.86 



Refinement quality 




x 2 

4.25 



Rp (%) 

3.99 



Rwp (%) 

5.48 




Lattice parameters 

EuCu 2 Sb 2 
This Work 

Ref. 39 

a (A) 

4.4876(2) 

4.504(1) 

c (A) 

10.7779(5) 

10.824(2) 

c/a 

2.4017(2) 

2.403(1) 

Well (A 3 ) 

217.05(3) 

219.58(14) 

Vm (cm 3 /mol) 

65.35 


Refinement quality 

x 2 

4.21 


R P (%) 

4.76 


7?wp (%) 

6.62 



TABLE II: Atomic coordinates and occupancies (Occ) ob¬ 
tained from the Rietveld refinements of powder XRD data 
for crushed crystals of EuCu 2 As 2 and EuCu 2 Sb 2 . Also shown 
are single-crystal data for the As z parameter from Ref. 39. 
The value for this parameter determined in the present work 
is denoted by “PW”. 


Atom 

Wyckoff 

symbol 

Occ 

(%) 

X 

y 

z 

(PW) 

z 

(Ref. 39) 

EuCu 2 As 2 

Eu 

2 a 

100 

0 

0 

0 

0 

Cu 

4 d 

98(1) 

0 

1/2 

1/4 

1/4 

As 

4e 

100 

0 

0 

0.3798(2) 

0.3762(1) 

EuCu 2 Sb 2 

Eu 

2 c 

100 

1/4 

1/4 

0.2381(3) 

0.2381(1) 

Cul 

2 a 

101 (1) 

3/4 

1/4 

0 

0 

Cu2 

2 c 

82(1) 

1/4 

1/4 

0.6366(9) 

0.6365(2) 

Sbl 

2b 

100 

3/4 

1/4 

1/2 

1/2 

Sb2 

2 c 

100 

1/4 

1/4 

0.8659(4) 

0.8692(1) 


each other, in contrast to only one type of stacked square 
lattice of Cu atoms in EuCu 2 As 2 . 

The crystallographic and refinement quality parame¬ 
ters obtained from the Rietveld refinements of the XRD 
data are listed in Tables I and II. While refining the XRD 
data the occupancies of Eu and As/Sb positions were 
kept fixed at their stoichiomentric values and the occu- 
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pation of the Cu site(s) were allowed to vary. The ther¬ 
mal parameters B were also fixed to B = 0 as there was 
no change in the lattice parameters and/or z-parameters 
within the error bars upon varying B. The occupancies 
are found to be stoichiomentric except for the Cu2 posi¬ 
tion in EuCu 2 Sb 2 which showed 18(1)% vacancies, for an 
overall composition EuCu 1 . 82 (i)Sb 2 . These occupancies 
are in agreement with our EDS results in Sec. II which 
suggested the presence of Cu vacancies in EuCu 2 Sb 2 but 
not in EUCU 2 AS 2 . Also shown in Tables I and II for com¬ 
parison are the lattice parameters and atomic positions 
for single crystals of both compounds reported in Ref. 39 
and the lattice parameters for a polycrystalline sample 
of EUCU 2 AS 2 reported in Ref. 40. 

The differences in the lattice parameters and unit 
cell volumes in Table I reported in the three studies of 
EuCu 2 As 2 and EuCii 2 Sb 2 are outside the respective er¬ 
ror bars stated, suggesting that the different samples of 
EUCU 2 AS 2 and EuCu 2 Sb 2 can have different composi¬ 
tions. This variability evidently arises from variations in 
the Cu site occupancies. Even though the unit cell vol¬ 
umes of EuCu 2 As2 and EuCu 2 Sb 2 differ by almost 20%, 
the c/a ratios are nearly identical at ~ 2.40. 


IV. PHYSICAL PROPERTIES OF EuCui 82 Sb 2 
CRYSTALS 

We begin with the physical properties of EuCui. 82 Sb 2 
because the interpretation of the magnetic data for this 
compound is more straightforward than for EUCU 2 AS 2 . 
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A. Magnetization and Magnetic Susceptibility 

The zero-field-cooled (ZFC) and field-cooled (FC) 
magnetic susceptibilities x = M/H of an EuCui. 82 Sb 2 
single crystal as a function of temperature T measured 
at different H aligned along the c axis (x c , H || c) and in 
the ab plane ( \ab, H T c) are shown in Figs. 3 and 4. At 
low H (e.g., at H = 0.01 T, Fig. 3) sharp anomalies are 
observed at T = 5.1 K in both \ab and Xc that we identify 
as the AFM ordering (Neel) temperature Tn. No hystere¬ 
sis is observed in ZFC and FC x(^) data [Fig. 3(b)]. An 
increase in H results in a shift of Tn to lower T (upper 
insets of Fig. 4), consistent with expectation for an AFM 
transition. The x(T < Tn) in Fig. 3 is anisotropic with 
Xab < Xc; indicating that the easy axis lies in the ab plane 
(collinear AFM ordering) or the easy plane is the ab plane 
(planar noncollinear AFM ordering). This easy-plane 
behavior is common for AFM ordering in the iron ar¬ 
senide family, e.g., in AFe 2 As 2 (A = Ca, Sr, Ba, Eu). 8,21 
As discussed in Sec. IV B below, the data in Fig. 3 are 
consistent with either a collinear A-type AFM structure 
with domains with orthogonal easy axes, or a planar non¬ 
collinear helix with the helix axis being the c axis. The 
A-type AFM candidate structure of EuCui. 82 Sb 2 for a 
single domain is illustrated in Fig. 2(b). 


FIG. 3: (Color online) Zero-field-cooled and field-cooled mag¬ 
netic susceptibility \ of EuCui. 82 Sb 2 as a function of temper¬ 
ature T in the temperature range 1.8-20 K measured on two 
different single crystals (crystal 1 and crystal 2) in a magnetic 
field H = 0.01 T applied along the c axis (xc, H |j c) and in 
the ab plane ( \ab , H _L c). 


In the paramagnetic (PM) state the x(T > Tn) data 
follow the Curie-Weiss behavior in Eq. (1). The lin¬ 
ear fits of x _1 (T) data (shown by straight red lines in 
Fig. 4) measured at H = 3 T in the temperature range 
50 K < T < 350 K by the Curie-Weiss law yield the Curie 
constants C and Weiss temperatures 9 p listed in Ta¬ 
ble III, along with the values of the parameter / = 0 p /Tn- 
Also listed are the effective magnetic moments calculated 
from the values of C using the relation p, e ff = V8C. 
The values of /x e fi are similar to the theoretical value 
/x e ff = g^/S(S + 1) hb = 7.94 g ,b for a free Eu +2 ion with 
spin S = 7/2 and spectroscopic splitting factor g = 2. 
This indicates that the Eu ions in EuCui. 82 Sb 2 are diva¬ 
lent. 

To further clarify the nature of the AFM structure, 
M(H) isotherms were obtained at eight temperatures 
between 1.8 and 300 K for H applied along the c axis 
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FIG. 4: (Color online) Zero-field-cooled inverse magnetic sus¬ 
ceptibility y _1 of an EuCui.8 2 Sb 2 single crystal (crystal 1) as 
a function of temperature T in the temperature range 1.8- 
350 K measured in a magnetic field of 3.0 T applied (a) in the 
afe-plane (Xab,H _L c) and, (b) along the c-axis ( Xc,H || c). 
The solid straight red lines are fits of the y _1 (T) data by 
the Curie-Weiss law (1) in the T range 50 K < T < 350 K. 
The upper insets in (a) and (b) show x(T) at low T and at 
different applied fields and the lower insets show fits of the 
Curie-Weiss law (1) to X _1 (T) data between 20 and 30 K in 
H = 0.1 T. 


(M c , H || c) and in the ab plane (M a b, H _L c) as shown in 
Figs. 5 and 6. Figure 5 shows that at 1.8 K, M increases 
almost linearly with H up to H ~ 3 T (except for a slope 
change near H = 0.5 T for H _L c), above which M tends 
to saturate with saturation moments = 6.66 /zb and 
Msat = 6-77 /zb at H = 5.5 T for M a b and M c , respec¬ 
tively. These saturation moments are close to the ex¬ 
pected /t sa t = gSg, b/Eu = 7/is/Eu assuming g = 2 and 
S = 7/2. No hysteresis was observed between the increas¬ 
ing and decreasing cycles of H for the M(H) isotherms 
at 1.8 K (data not shown). At 5 K « Tn = 5.1 K, the 
M(H) curves in Fig. 5 for the two held directions exhibit 
very similar behaviors, as expected for the PM state. At 


TABLE III: Antiferromagnetic ordering temperature Tn and 
the parameters obtained from Curie-Weiss fits of the magnetic 
susceptibility data for EuCu 2 As 2 and EuCui.s 2 Sb 2 , where C 
is the Curie constant and 6 P is the Weiss temperature. The 
parameter / is defined as / = 0 p /Tn. The effective magnetic 
moment per Eu atom, /i e ff, is obtained from C according to 
Peff = y/8C. The H = 3 T data were fitted by the Curie- 
Weiss law from 50 to 350 K, whereas the H = 0.1 T data were 
fitted from 20 to 30 K to obtain more accurate estimates of 
the Weiss temperatures. 


Compound 

Field 

direction 

Tn 

(K) 

C 

/ cm 3 K \ 

V mol / 

Or 

(K) 

/ 

fte ff 

\ En ) 

H = 3 T 
EuCu 2 As 2 

H 1 c 

17.5 

7.44(1) 

+19.0(2) 

1.09 

7.72(1) 


H ||c 

17.5 

7.64(3) 

+17.2(5) 

0.98 

7.82(2) 

EuCui. 8 2Sb 2 

H T c 

5.1 

7.41(1) 

+2.2(4) 

0.43 

7.70(1) 


H || c 

5.1 

7.54(1) 

+2.1(1) 

0.41 

7.77(1) 

H = 0.1 T 
EuCui. 82 Sb 2 

H _L c 

5.1 

7.62(1) 

+ 1.26(4) 

0.25 

7.81(1) 


H\\c 

5.1 

7.71(1) 

+ 1.69(4) 

0.33 

7.85(1) 


temperatures T > 50 K, M is nearly proportional to H 
as expected in the PM regime with gg^H / (k-g,T) <C 1 for 
H < 5.5 T. 

We obtained high-resolution M(H) data at 1.8 K with 
H _L c for H < 0.75 T as shown in Fig. 6. One observes a 
nonlinear behavior of M(H) in this held range. To study 
this nonlinearity in detail, shown in the inset of Fig. 6 are 
plots of dM/dH versus H for H _L c at 1.8 and 5 K and 
for H || c at 1.8 K. For H || c at 1.8 K, dM/dH versus 
H is featureless, confirming that the ab plane is the easy 
plane and the c axis is a hard axis. The dM/dH versus 
H data at 5 K for H _L c is also featureless, as expected 
for the PM state. 

On the other hand, the dM/dH versus H data at 1.8 K 
for H _L c in the inset of Fig. 6 exhibit two important fea¬ 
tures. First x = dM/dH is constant for 0 < H < 0.1 T. 
Second, for H > 0.1 T, dM/dH increases and exhibits 
a pronounced peak at H « 0.5 T. These two features 
suggest a spin-hop transition that is distributed in held. 
In such a transition, the ordered moments hop from an 
orientation that is not perpendicular to H to a perpendic¬ 
ular orientation. The held scale over which the magneti¬ 
zation changes sharply in the inset of Fig. 6 (< 0.7 T) at 
T = 1.8 K and H T c is of the order expected from a spin- 
hop transition associated with magnetic dipole interac¬ 
tion anisotropy as deduced from the data in Table IV 
below, where the ordered moment direction hops from 
the low-susceptibility direction parallel to the ab plane 
to the high-susceptibility direction approximately paral¬ 
lel to the c axis (see Fig. 3). At higher helds the ordered 
moments increasingly point towards the applied held di¬ 
rection until they are parallel to the applied held with 
parallel moment approximately equal to the saturation 
moment of ~ 7 /zb/Eu, which happens at the critical 
held H c ~ 3 T as seen in Fig. 6, which is discussed in 
quantitative detail in the following section. For H > H c 














































H (T) 

FIG. 5: (Color online) Isothermal magnetization M of an 
EuCui. 82 Sb 2 single crystal (crystal 1) as a function of internal 
magnetic field H measured at the indicated temperatures for 
H (a) in the ab plane ( M a i,H _L c) and (b) along the c axis 
(M c ,H || c). 


FIG. 6: (Color online) Isothermal magnetization M of an 
EuCui. 8 2 Sb 2 single crystal (crystal 2) as a function of internal 
magnetic field H measured at 1.8 K for H J_ c. Inset: The 
derivative dM/dH versus H at 1.8 K for H _L c and H || c, 
and at 5 K for H _L c. 


bors. 

Here we discuss two candidates for the magnetic struc¬ 
ture. The first is multiple domains of a collinear A-type 
AFM with the ordered moments aligned in the ab plane, 
and the second is a planar noncollinear helical AFM 
structure with the ordered moments again aligned in the 
ab plane which is perpendicular to the helix c axis. A fig¬ 
ure showing the helical AFM structure is given in Ref. 23. 
Here we only consider Heisenberg exchange interactions 
between the spins. We discuss later how some of the pre¬ 
dictions are modified by the presence of magnetic dipole 
interactions. 


1. Magnetic Susceptibility: Collinear Antiferromagnetic 
Ordering 


the system is in the PM phase. 


B. Molecular Field Theory Analysis of Magnetic 
Properties 

Because the Eu spin S = 7/2 is large, the quantum 
fluctuations associated with finite spin should be small 
and we expect molecular field theory (MFT) to fit the 
experimental magnetic susceptibility and heat capacity 
data rather well. We use a version of MFT for an- 
tiferromagnets developed by one of us for systems of 
identical crystallographically-equivalent spins interacting 
by Heisenberg exchange that does not use the concept 
of magnetic sublattices. 22-24 Instead, the magnetic and 
thermal properties are calculated solely from the ex¬ 
change interactions of an arbitrary spin with its neigh- 


In the Weiss MFT for a system of identical crystal- 
lographically equivalent spins interacting by Heisenberg 
exchange with Hamiltonian fi = 'Yh <l]> • S 7 , where 

the sum is over distinct pairs of spins interacting with 
exchange constants Jij and a positive Jij corresponds to 
an AFM interaction and a negative one to a FM interac¬ 
tion, the magnetic susceptibility X|| parallel to the easy 
axis of a collinear AFM at T < Tn is given by the law of 
corresponding states 23,24 


where 


/ 


0P_ 

T n ’ 


XII ( T ) = 1-/ 

xCTn) T*(i)-/’ 


(5a) 


t 


T 


r*{t) 


{S + l)t 
3 B' s {y Q y 


3^o 

(5 + 1)*’ 
(5b) 


















the ordered moment versus temperature in zero field 
is denoted by /j, 0 , the reduced ordered moment JIq = 
Mo/Vsat = Mo/(<?5^b) is determined by solving 

Mo = B s {yo), (5c) 


B' s (yo ) = [dBs(y)/dy\\y — y0 , and our unconventional def¬ 
inition of the Brillouin function Bs(y) is 


Bs(y) 


1 

25 


| (25 + 1) coth (25 + 1) 


2 J 



(5d) 

Within MFT, the susceptibility in the ordered state of 
an AFM with the field applied perpendicular to the easy 
axis of a collinear AFM or to the easy plane of a planar 
noncollinear AFM is independent of T, i.e., 


X±(T 1 <T N ) = X (T N ). (6) 


In the case of A-type AFM ordering with the ordered 
moments aligned in the ab plane, there are two equiv¬ 
alent orthogonal directions for the easy axis due to the 
two equivalent orthogonal a and b axes of the tetragonal 
unit cell. This gives rise to two equivalent A-type AFM 
domains with easy axes in the ab plane that are orthogo¬ 
nal to each other. The fractional populations of the two 
domains need not be the same. Let x be the fractional 
population of domains with the easy axis perpendicular 
to the applied field with H + c. Then the fractional 
population of domains with the easy axis parallel to the 
applied field is 1 —x. The average (measured) susceptility 
in the ab plane is then 


Xab ave = +XX + (1 ~ x)x\\ (T) = X\(T N ) + (1 - x)x\\ (T), 

(7a) 


where we used Eq. (6) for x_l- Since xii^ = 0) = 0, 
the value of x is given by 

_ Xabave(T = 0) 

r “ x(Tn) ■ 


23,24 

(7b) 


From the data in Fig. 7 one obtains x = 0.56. The 
additional parameters needed to fit the Xab(T ) data for 
T < Tn by the MFT are /, the saturation moment fj , SSLt = 
gSy, b and x(2n)- For EuCui.g 2 Sb 2 and H + c we use / = 
0.25 from Table III, /r sa t = gSy b = 7 Pb using g = 2 and 
5 = 7/2, and x(Fn) = 1-04 cm 3 /mol from Fig. 7. The fit 
of the experimental Xab{T) data for T < Tn by Eqs. (5a) 
and (7) is shown by the solid blue curve in Fig. 7(a) with 
no adjustable parameters. The experimental data are 
seen to be well represented by the MFT prediction for 
collinear A-type AFM ordering with, on average, 56% of 
the AFM domains having the easy axis in the ab plane 
oriented perpendicular to the applied field in the ab plane 
and 44% of them having the easy axis parallel to the 
applied field. 


2. Magnetic Susceptibility: Planar Noncollinear Helical 
Antiferromagnetic Ordering 

The in-plane susceptibility Xxy for a planar non¬ 
collinear helical AFM system is given within MFT by 23,24 




FIG. 7: (Color online) Fits of the Xab (T < Tn) data with 
H + c (crystal 2) shown in Fig. 3 for EuCui. 82 Sb 2 by MFT 
for (a) a multi-domain A-type AFM structure and for (b) a 
c-axis helical AFM structure. The fits are about equally good 
and hence the fits do not distinguish between the two AFM 
structures. 


Xxy(T < Tn) _ (1 + r* + 2/ + 4£?*)(1 — /)/2 

X(Tn) ~ (r* + B*)( 1 + B *) - (/ + B*Y ’ 
where 

B* = 2(1 — /) cos [kd) [1 + cos {kd)\ — f (8b) 

and kd is turn angle in radians between the ordered mo¬ 
ments in adjacent layers along the helix axis. For T = 0 
one obtains the simple result 

Xxy{T = 0) = _1_ 

x(7n) 2[l + 2 cos (kd) + 2 cos 2 (fed)] ’ 

which only depends on kd. For 0.5 < Xxy{T = 
0)/x(7n) < 1 there are two solutions for kd in 

Eq. (9). 23,24 
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Taking Xxy(T = 0 )/x(Tn) = 0.56 from Fig. 7 and solv¬ 
ing Eq. (9) for kd gives the two solutions kd = 93° and 
kd = 161° for the turn angle of the helix between ad¬ 
jacent FM-aligned layers along the helix axis in the or¬ 
dered state of EuCui. 82 Sb 2 at T = 0. Both solutions 
correspond to dominant AFM interactions between an 
ordered moment in one layer and the ordered moments 
in either of the two adjacent layers because the projec¬ 
tion of an ordered moment in one layer on an ordered 
moment in an adjacent layer is negative. The predic¬ 
tions of Xx V {T < T n )/ X (T n ) for the two values of kd are 
the same and therefore no decision as to which angle is 
more appropriate is possible from fitting the x(T) data 
alone. The planar noncollinear helical x(T) according to 
Eqs. (8) for S = 7/2 and / = 0.25 with a turn angle 
kd = 93° or 161° is shown in Fig. 7(b). The x(T < In) 
data are described by the helical AFM model and the 
A-type AFM model equally well. 

The magnetic structure for this helical model can be 
visualized from Fig. 2(b) with the difference that now the 
angle between the moments of adjacent layers along the 
c axis is 93° or 161° as opposed to 180° for the A-type 
AFM model. However, from a comparison of Figs. 7(a) 
and 7(b), one sees that when A-type AFM domains with 
orthogonal easy axes within the ab plane are present in 
a compound (which make an angle of 90° to each other 
in the two types of domains), the temperature depen¬ 
dence of the in-plane Xxy is very similar to that of the 
corresponding helical AFM with a turn angle of 93° (or 
161°). 


3. Influence of Magnetic Dipole Interactions 


The response of a magnetic moment in a sample is 
determined by the local effective magnetic induction it 
sees. After crystal shape effects are accounted for, as is 
consistently done in the M(H) and x(T) data presented 
throughout this paper, one has 

A -TT 

B\Tai = H Qa + Y M a + ( 10 ) 

where the second term is the macroscopic Lorentz field 
inside a spherical Lorentz cavity, M a = n/V s pin is the 
magnetic moment per unit volume, and the volume per 
spin is E S pi n = Fceii/2 = a 2 c/ 2. The third term is 

the contribution to the local magnetic induction in the 
a th direction due to the discrete point magnetic dipoles 
around a central spin within the Lorentz cavity centered 
on the central spin given by 47,48 


B 


near 
int a i 




(ii) 


where the factor of two is necessary in the first equal¬ 
ity because Ei is evenly split between the central mo¬ 
ment and a neighbor, whereas Bf/f a / i arises only from 
the neighbor and Aka is the eigenvalue of the magnetic 


dipole interaction tensor Gj(k) defined as follows. The 
eigenenergies Ei are 

Ei = -e /t?Gi(k)/ij, (12a) 


where 


e = 



(12b) 


has dimensions of energy, a is the tetragonal a-axis lat¬ 
tice parameter of a simple tetragonal or body-centered 
tetragonal spin lettice and 


G*(k) = £7- 




(rji/a) 5 


e 




(12c) 


is a dimensionless interaction tensor for collinear mag¬ 
netic ordering, r$ is the position in Cartesian coordi¬ 
nates of a central spin i at which the net magnetic in¬ 
duction due to spins at positions r 7 is calculated, r ji = 
Yj — r*, rji = |rjj| and k is the magnetic wave vector. La¬ 
beling the eigenvalues of Gb(k) as Aka, Eq. (12a) gives 
the eigenenergies 


Eia — —e Aka, (lid) 

and the eigenvectors fi are the ordered moment axes of 
the collinear magnetic structure. For FM alignment, 
which can include either an ordered FM structure or an 
induced alignment due to an applied magnetic field, one 
has k = (0,0,0) and for AFM wave vectors the com¬ 
ponents of k are expressed in conventional primitive- 
tetragonal reciprocal lattice units (rlu) 2ir/a and 2tt/c. 

Shown in Table IV are the eigenvalues Aka and eigen¬ 
vectors fi of Gj(k) calculated by direct lattice summation 
for a variety of moment configurations for body-centered- 
tetragonal spin lattices with c/a = 2.40 as found for the 
Eu positions in EUCU 2 AS 2 , EuCui. 82 Sb 2 and for the Gd 
positions in GdAu 2 Si 2 - 48 For a moment /i = gS[iB with 
g = 2 and S = 7/2 for Eu +2 or Gd +2 and the a-axis lat¬ 
tice parameters for EUCU 2 AS 2 and EuCui. 82 Sb 2 in Ta¬ 
ble I, Eq. (12b) gives 

e = 17.34 /reV = 0.2012 K (EuCu 2 As 2 ), (13a) 
e = 14.55 /reV = 0.1688 K (EuCui. 82 Sb 2 ). (13b) 

Using these values, and assuming stoichiomentric formu¬ 
las EuCu 2 As 2 and EuCu 2 Sb 2 , the values of E ia for each 
of the k values shown were calculated from Eq. (12a) and 
are listed in Table IV in units of both /xeV and K. The 
corresponding values of B/ aa D obtained using Eq. (11) 
are also shown for the lowest-energy eigenvector for 
each k. Our eigenvectors for k = (A,0,A) agree with 
those calculated in Ref. 48, but our Ei a values are about 
a factor of 5-6 larger in magnitude than in Ref. 48 for 
unknown reasons. 

One sees from Table IV that the magnetic dipole inter¬ 
action favors either k = (1,0, 0) rlu or k = (1,0, rlu, 
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TABLE IV: Eigenvalues Aka and eigenvectors (i = [g x , g y , g z \ of the magnetic dipole interaction tensor Gi(k) in Eq. (12c) 
for c/a = 2.40 and six values of the magnetic wave vector k for collinear magnetic order. The accuracy of the Aka values is 
estimated to be ±0.001. The Ei values for EUCU 2 AS 2 and EuCui.8 2 Sb 2 were obtained from Eq. (12d) using the Aka values 
in the second column and the e values in Eqs. (13) which assume that the ordered moment magnitude is g = 7/xb/Eu. Also 
shown are local magnetic induction values B/// 1 / = —2Ei/(7g,B) obtained from Eq. (11). 


k 

Aka 

A 

-Ei 

EuCu 2 As 2 

(peV) 

—Ei/kB 

EuCu 2 As2 

(K) 

pnear 
-*-'int a i 

EuCuzAsa 

(G) 

-Ei 

EuCui.82Sb2 

(peV) 

—Ei/kB 

EuCui.8 2 Sb 2 

(K) 

pnear 
-^int a i 

EuCui.8 2 Sb 2 

(G) 

o, 0) 

5.100 

[0, 1 , 0] 

88.5 

1.026 

4366 

74.2 

0.861 

3664 


0.977 

[0, 0, 1] 

17.0 

0.197 


14.2 

0.165 



-6.077 

[1, 0, 0] 

-105.4 

-1.223 


-88.5 

-1.026 


(§, 0, i) 

5.100 

[0, 1 , 0] 

88.4 

1.026 

4364 

74.2 

0.861 

3664 


1.328 

[0.2366,0,-0.9716] 

23.0 

0.267 


19.3 

0.224 



-6.428 

[-0.9716,0,-0.2366] 

-111.5 

-1.293 


-93.6 

-1.085 


(0, 0, i) 

4.517 

[1, 0, 0] 

78.3 

0.909 

3866 

65.7 

0.763 

3244 


4.517 

[0, 1 , 0] 

78.3 

0.909 


65.7 

0.763 



-9.035 

[0, 0, 1] 

-156.7 

-1.818 


-131.5 

-1.525 


(0, 0, 1) 

4.439 

[1, 0, 0] 

77.0 

0.893 

3798 

64.6 

0.749 

3188 

(A-type AFM} 

4.439 

[0, 1 , 0] 

77.0 

0.893 


64.6 

0.749 



-8.877 

[0, 0, 1] 

-153.9 

-1.786 


-129.2 

-1.499 


(i i °) 

2.652 

[0, 0, 1] 

46.0 

0.534 

2270 

38.6 

0.448 

1906 


-0.789 

[72’”72’ 0 ] 

-13.7 

-0.159 


-11.5 

-0.133 



-1.863 

[ 72 ’ 75’°J 

-32.3 

-0.375 


-27.1 

-0.315 


(§, i 1) 

2.640 

[0, 0, 1 ] 

45.8 

0.5311 

2260 

38.4 

0.446 

1896 

(G-type or 

-1.320 

[1,0,0] 

-22.9 

-0.266 


-19.2 

-0.223 


Neel-type AFM} 

-1.320 

[0,1,0] 

-22.9 

-0.266 


-19.2 

-0.223 


(0, 0, 0 ) 

1.105 

[1, 0, 0 ] 

19.2 

0.222 

946 

16.1 

0.187 

794 

(ferromagnetic 

1.105 

[0, 1, 0] 

19.2 

0.222 


16.1 

0.187 


alignment} 

-2.210 

[0, 0, 1 ] 

-38.3 

-0.445 


-32.2 

-0.373 



neither of which corresponds to the AFM structures pos¬ 
tulated above for EuCui. 82 Sb 2 . However, the ordered 
moment direction in both cases is along the b axis, in 
agreement with our measurements to be in the ab plane. 
The FM structure has the highest energy among the k 
values listed in the table. 

For the collinear A-type AFM structure possibibility 
with k = (0,0,1) rlu discussed in Sec. IV B 1, the lowest- 
energy eigenvector for this k value is also in the ab plane, 
consistent with the susceptibility data and fit in Figs. 3 
and 7, respectively. On the other hand, this k does 
not have the lowest ordering energy among the k values 
listed, which means that Heisenberg exchange interac¬ 
tions instead of dipole interactions determine the AFM 
structure, but with the dipole interaction determining or 
at least contributing to the easy axis for the ordering. Of 
our two postulated AFM structures, neutron diffraction 
measurements have determined that EuCu 2 Sb 2 orders in 
the A-type AFM structure with k = (0,0,1) and the easy 
axis indeed lies in the ab plane as predicted for the mag¬ 
netic dipole interaction. 49 

Using MFT, the Neel temperature Tna for A-type 
AFM ordering and g = (1,0,0) arising from the 
anisotropic magnetic dipole interaction is 22-24 


Setting g = 2, S = 7/2, A ka = 4.439 from Table IV and 
a = 4.4876 A for EuCui .82 Sb 2 from Table I, one obtains 

T na = 0.64 K (EuCui. 82 Sb 2 ). (15) 

Because the local fields due to the exchange and dipo¬ 
lar interactions are additive in their contributions to Tn 
within MFT, this is the amount by which Tn increases 
due to the magnetic dipole interaction. Since the mea¬ 
surements give Tn = 5.1 K, this is only a 13% effect, 
which indicates that the measured Tn is mainly due to 
exchange interactions. 

In the PM state above the Neel temperature, all mo¬ 
ments are aligned in the same direction a as the applied 
magnetic field H$ a [the magnetic wave vector is k = 
(0,0,0)]. Then in the limit of low field and after cor¬ 
rection for shape effects leading to a demagnetizing field, 
MFT for dipolar interactions yields the Curie-Weiss law 


where the single-spin Curie constant 50 C\ and the Weiss 
temperature 9 P a for a bet spin lattice are 


Tna 


g 2 S(S + l)^|A ka 
3 a 3 ks 


(14) 


, = g 2 s{s + 1)/4 

3A;b 


(17a) 
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0pA ~&(3^ + A(o ’°’ o) “)- (1?b) 

Using g = 2, S = 7/2, a = 4.4876 A and c/a = 2.40 
from Table I, and A( 0 ,o,o)[i,o,o] = 1-105 and A( 0 ,o,o)[o,o,i] = 
—2.210 from Table IV, Eqs. (17) give 

d pA = 0.59 K (H ic), (18a) 

d pA = 0.35 K (H || c). (18b) 

The contributions of the magnetic dipole interaction to 
d pA for H 1 c and H || c are both positive and hence 
ferromagneticlikc. 


4. Exchange Constants and Exchange Field 

One can estimate the exchange interactions within 
MFT. In this discussion we correct for the contributions 
of magnetic dipole interactions to the observed values of 
Tn and 6 P . According to Eqs. (15) and (18b), the Tn is 
increased by 0.64 K, 9 p (H _L c) is increased by 0.59 K 
and 9 p (H || c) is increased by 0.35 K due to the mag¬ 
netic dipole interactions. Subtracting these values from 
the observed Tn = 5.1 K, 9 p (H J_ c) = 1.26 K and 
9 p (H || c) = 1.70 K in Table III gives the contributions 
due to the exchange interactions Jij as 


T nj = 4.5 K, (19a) 

Q p j(H 1 c) = 0.67 K, 6 pJ {H || c) = 1.35 IC. (19b) 

Thus there appears to be another source of anisotropy 
present in addition to the magnetic dipole interaction. 
Taking the spherical average of these two values as an 
approximation gives 


6 p j = 0.90 K. 


( 20 ) 


The value fj arising from the exchange interactions as 

fj=?r~= 0 - 20 . ( 21 ) 

Jn j 

In MFT 9 p j and fj are related to the exchange inter¬ 
actions by 23,24 


Xnj 


S(S+ 1) 
3fc B 


/ ( Jij COS (f)j j , 
j 


( 22 ) 


and 


0 P J = ~ 


S{S+ 1) 
3 fc B 




fj 


E, Jij COS (f>ji 


(23) 


(24) 


where the sums are over all neighbors j with which a cen¬ 
tral spin i interacts with respective exchange constant Jij 
and (i>ji is the angle between ordered moments jij and /q 
in the magnetically-ordered state. For the A-type AFM 
structure in Fig. 2(b) with a fourfold in-plane nearest- 
neighbor interaction J\ and twofold interlayer interaction 
J c these expressions become 

Tnj = - 5( 1 +1) (4J 1 -8J c ), (25) 

oKb 


0 p J = ~ S ^t + ^ QJi+SJc) (26) 

3 k B 

and 


_ 9 P j _ Ji+2J C 
Ji - 2 J c 


(27) 


Using fj = 0.20, 6 p j = 0.90 K and S = 7/2, from 
Eqs. (26) and (27) we obtain 


^-= -0.129 K, —^ = 0.043 K. (28) 
kb kb 

Thus J\ is FM and J c is AFM, consistent with the A-type 
AFM structure in Fig. 2(b). 

The exchange field -ff eX cho can be estimated from the 
values of J\ and J c using the relation 23,24 


-ffexch ? 


2 .,2 J j j Mi cos 0i i ? 

9 Mb j 


which for the present case is 

H exch0 (T = 0) = —^j-[4Ji - 2 J c \, 
M 2 Mb 


(29) 


(30) 


where g = 2 and go = 7 fiB- This gives -ffexcho^ = 
0) = 22.4 kOe. This is about a factor of 7 larger than 
the dipolar magnetic induction of 3188 G for k = (0,0, 1) 
in Table IV, as expected since the exchange interactions 
mainly determine Tn as noted above. 


5. High Magnetic Fields Perpendicular to the Easy Axis or 
Plane of the Antiferromagnetic Structure 

The version of MFT considered in this paper 23,24 does 
not utilize the concept of magnetic sublattices, so the 
thermal and magnetic behavior of each ordered moment 
in H = 0 is the same and only depends on its interactions 
with its neighbors. Because the response of each ordered 
moment to a field applied perpendicular to the easy axis 
of a collinear AFM or to the easy plane of a planar non- 
collinear AFM is the same since the former structure is a 
special case of the latter, the same MFT can be applied 
to obtain a law of corresponding states for a given S for 
the response of a given moment to an applied perpendic¬ 
ular field for both types of structures. We will use these 
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FIG. 8: (Color online) Magnetization M for crystal 2 of 
EuCui.gzSba, normalized by its saturation value M a = 
NgSg B , versus internal magnetic field H for H || c for sev¬ 
eral temperatures T < Tn = 5.1 K. This field direction is 
perpendicular to the easy axes or plane of the two candidate 
AFM structures considered in this paper. Also shown as solid 
curves are fits of the data by MFT at the respective temper¬ 
atures according to Eqs. (34). 


MFT results 24 to describe the high-field perpendicular 
magnetization of EuCui. 82 Sb 2 at T < Tw¬ 
in MFT, the initial slope of the magnetic moment per 
spin fi± of a collinear or planar noncollinear AFM versus 
perpendicular magnetic field H is independent of T for 
0 < T < Tn and the M(H) is given by 

Mx = X_l H, x± = x( t n), (31) 


where x(Tn) is the low-field single-spin susceptibility at 
T = Tn. The perpendicular critical field H c ± is the 
perpendicular field at which a collinear or planar non- 
collincar AFM undergoes a second-order transition from 
the (canted) AFM state to the paramagnetic (PM) state. 
As long as H < H c ±, the ordered moment magnitude is 
independent of field even though the ordered moments 
are progressively tilting towards the field with increasing 
field, and the reduced ordered moment is therefore given 
in this field range by the expression for the zero-field re¬ 
duced ordered moment flo in Eq. (5c). The critical field 
is the field at which the angle between the ordered mo¬ 
ments and the perpendicular field becomes equal to zero 
and hence Eq. (31) gives 


H c± (T = 0) 


N A fi 0 

x(r N ) ’ 


(32) 


where N A is Avogadro’s number, = gSg, B and x(Tn) 
is here expressed per mole of spins. For H > H c ± the sys¬ 
tem is in the PM state and hence the magnetic moments 
in the direction of the field are considered to be purely 
field-induced moments, as opposed to ordered moments 
present in the AFM state. Taking g = 2 and S = 7/2 ap¬ 
propriate to Eu +2 spins and %(Tn) = 1.3 cm 3 /(mol Eu) 


from Fig. 3(b) gives H C ±(T = 0) ~ 26 kOe, similar to 
the value of ~ 30 kOe obtained from Fig. 5. 

It is convenient to define the reduced parameters 


= gg, B H 
k B T N 


h c ± 


gH B H c _ l 

fc B T N 


Mx = —, (33) 

M sat 


where /u sa t = gSfi B . Then the H dependence of the or¬ 
dered plus induced moment per spin /zj_ for fields applied 
perpendicular to the easy axis (collinear AFM) or easy 
plane (planar noncollinear AFM) of a Heisenberg AFM 
is described by the same law of corresponding states for 
a given S, given by 24 


MX 


5 + 1 
3(1-/) 


h 


Mx = B s 


' 3/mx 
,(5+l)i 


h 

t 


(h < h c ±), 

(34) 

(h X h c _ i_), 


where fix = M/M s , M = M s = N/j, sat , N is the 

number of spins in the system and t = T /Tn. 

The normalized M(H)/M S data for EuCui.s 2 Sb 2 crys¬ 
tal 2 at T = 1.8, 3.0 and 4.0 K are shown for H \\ c 
by the open symbols in Fig. 8. Also shown are the pre¬ 
dictions of MFT for M(H)/M b using Eqs. (34) and the 
value / = 0.4 for H \\ c with H = 3 T from Table III. 
The theory is seen to be in semiquantitative agreement 
with the experimental data. 


C. Heat Capacity 

The heat capacity at constant pressure C p for an 
EuCui. 82 Sb 2 crystal as a function of T is shown in Fig. 9. 
A pronounced A-type anomaly at 5.1 K (inset of Fig. 9) 
is observed that confirms the intrinsic nature of AFM 
ordering in this compound. Due to the strong influence 
of the AFM order below Tn and short-range AFM order 
above Tn it is difficult to estimate the electronic contri¬ 
bution to C P (T). The strong magnetic contribution at 
low temperatures prevented us from carrying out a con¬ 
ventional C/T versus T 2 fit to the C/T versus T 2 data 
to obtain the Sommerfeld coefficient 7 and the Debye T 3 
coefficient /?. 

From Fig. 9, at T = 300 K the C p attains a value 
of ss 123 J/molK which is close to the expected clas¬ 
sical Dulong-Petit value Cy = 3 nR = 14.461? = 
120.2 J/molK at constant volume , 50,51 where n = 4.82 
is the number of atoms per formula unit (f.u.) and R is 
the molar gas constant. We analyzed the C P (T) data in 
the PM state within the framework of the Debye lattice 
heat capacity model. We fitted the C P (T) data by 

C p (T) = 7 T + nCy Debye (T), (35a) 

where C\ Debye (T) represents the Debye lattice heat ca¬ 
pacity. The Debye model describes the lattice heat ca¬ 
pacity due to acoustic phonons at constant volume V 
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FIG. 9: (Color online) Heat capacity C p of an EuCui. 82 Sb 2 
single crystal as a function of temperature T in the temper¬ 
ature range 1.8-300 K measured in zero magnetic field. The 
solid curve is the fitted sum of the contributions from the 
Debye lattice heat capacity Cy Debye (T) and electronic heat 
capacity 7 T according to Eq. (35a). Inset: Expanded view of 
the low-T C P (T) of EuCui. 82 Sb 2 , of the nonmagnetic refer¬ 
ence compound SrCu 2 Sb 2 , 38 and of estimated lattice contri¬ 
bution after correcting for the difference in formula masses of 
EuCu 2 Sb 2 and EuCui .82 Sb2. 


which is given per mole of atoms by 51 

/ T \ 3 / ,e>D / T x 4 ’e x 

Cyr>ebye{T) = 9R J ( e x _ 1)2 ' (35b) 

The solid curve in Fig. 9 represents the fit of the C p (T ) 
data for 15 K < T < 300 K by Eqs. (35) which is ob¬ 
tained using the analytic Pade approximant fitting func¬ 
tion for Cy Debye (P) • 52 The fit gave the Sommerfeld coef¬ 
ficient 7 = 17(2) mJ/molK 2 and the Debye temperature 
@ D = 196(2) K. 

Now we estimate the magnetic contribution C' mag (T) 
to the heat capacity of EuCui. 82 Sb 2 . For this we used the 
C P (T) data of SrCu 2 Sb 2 to estimate the lattice contribu¬ 
tion to the heat capacity of EuCui. 82 Sb 2 . SrCu 2 Sb 2 also 
forms in the same CaBe 2 Ge 2 -type primitive tetragonal 
structure as EuCui. 82 Sb 2 . 38 However, they have different 
formula masses, therefore the C P (T ) data of SrCu 2 Sb 2 
need to be corrected for the mass difference in order to 
obtain an accurate estimate of the lattice heat capacity 
of EuCui. 82 Sb 2 . As seen from Eq. (35b), the lattice heat 
capacity is a function of T/0 d and 0 d depends on the 
formula mass M according to ©d ~ 1/M 1 / 2 , or equiv¬ 
alently, T/0 d ~ M 1 / 2 . Therefore the measured T of 
SrCu 2 Sb 2 was scaled as 

T 

rj p* /gg\ 

(M E uCu 2 Sb2 / A/srCu 2 Sb 2 ) ^ 

The mass-corrected lattice contribution thus obtained for 
EuCui. 82 Sb 2 is shown in the inset of Fig. 9. The C' mag (T) 



FIG. 10: (Color online) (a) Magnetic contribution Cmag to the 
heat capacity of EuCui. 82 Sb 2 plotted as C mag (T)/T versus T. 
The solid curve represents the MFT prediction for S = 7/2 
and Tn =5.1 K. (b) Magnetic contribution S mag to the en¬ 
tropy versus T. Shown are the data before (red empty circles) 
and after (black empty squares) formula-mass correction for 
the lattice heat capacity. 


of EuCui. 82 Sb 2 obtained by subtracting the lattice con¬ 
tribution from the measured C P (T) data of EuCui. 82 Sb 2 
is shown in Fig. 10(a) as C' mag (T)/T versus T, where 
a sharp anomaly due to the AFM transition is evident. 
Further we observe that the C' mag (T) is nonzero even 
above the Tn and becomes negligible above about 15 K. 
This nonzero contribution to C mag for Tn < T < 15 K 
reflects the presence of dynamic short-range AFM corre¬ 
lations above Tn- 

In an AFM state, spin-wave theory predicts a T 3 be¬ 
havior of C' mag (T) for a three-dimensional (3D) AFM 
with negligible anisotropy gap or a T 2 dependence in 
(quasi-) 2D. However, below Tn there is significant 
deviation from the expected power-law behavior. In¬ 
stead, a plateau is observed in C' mag (T)/T versus T 
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[see Fig. 10(a)]. This plateau in C mag (T)/T or a cor¬ 
responding broad hump in C' mag (T) below T n has been 
observed in several Eu +2 and Gd +3 compounds, such 
as in Gd 2 Fe 3 Si 5 (Ref. 53), GdCu 2 Si 2 (Ref. 54), EuB 6 
(Ref. 55), EuRh 2 Si 2 (Ref. 56), and EuCo 2 Ge 2 (Ref. 57). 
Attempts have been made in past to explain the origin 
of the hump in C P (T) below Tn. 53,55,58,59 j n orc j e r to 
explain the T-linear behavior of ordered state C P (T)/T 
data of Gd 2 Fe 3 Si 5 Vining and Shelton 53 speculated that 
this kind of behavior may be due to low-dimensional spin- 
wave type ordering. Siillow et al. 55 interpreted the obser¬ 
vation of plateau in C P (T)/T of EuBg in terms of split¬ 
ting of ground state multiplet of Eu +2 by internal mag¬ 
netic field. According to Fishman and Liu, 58 within the 
Heisenberg model quantum fluctuations can give rise to 
a hump in C P (T) of a FM system below Tn. 

A hump in C mag (T) arises naturally in the MFT of the 
ordered-state heat capacity of system having a (2 S + 1)- 
fold degenerate ground state for large S'. 22 ’ 59 The origin 
of the hump can be understood in terms of the entropy 
that increases with increasing S. In order to accommo¬ 
date the increased entropy a hump appears in the heat 
capacity. The temperature T* at which hump devel¬ 
ops decreases with increasing S. As shown in Ref. 22, 
for S = 7/2 there appears a hump in C mag (T) data 
at T* < Tn/3. The solid curve in Fig. 10(a) repre¬ 
sents the mean-field theoretical C mag (T)/T calculated for 
Tn = 5.1 K (Ref. 22) which reproduces the trend of the 
data. The reduced experimental value of C' mag (T)/T be¬ 
low Tn is due to the presence of short-range ordering 
above Tn as discussed next. 

The experimental magnetic contribution to the entropy 
Smag (T) was estimated by integrating the C' mag (T)/T 
versus T data according to 

S mag (T) = £ Cma z {T) dT. (37) 

The T dependence of S mag for EuCui. 82 Sb 2 is shown in 
Fig. 10(b) for the temperature range 0 K < T < 30 K. 
The S mag (T ) between 0 and 1.8 K was obtained by ex¬ 
trapolating C' mag (T)/T to T = 0 assuming the MFT be¬ 
havior as shown by the dotted curve in Fig. 10(a). It is 
seen from Fig. 10(b) that if the S mag (T ) is uncorrected 
for the lattice contribution, the obtained .Sjnag exceeds 
the expected Sjnag = .Rings' + 1) = Tin8 for a divalent 
Eu (S = 7/2). This illustrates that the lattice contribu¬ 
tion must be corrected for to obtain accurate estimates 
of the magnetic entropy above a few Kelvins. The cor¬ 
rected S mag attains a value of 16.5 J/mol K at 10 K which 
is 95% of the expected high-T limit 7? In 8. The mag¬ 
netic entropy of Rln8 is fully recovered at 20 K. Thus, 
consistent with the x(T) data, the S mag (T) data show 
that the Eu atoms are in the Eu +2 state with S = 7/2 in 
EuCui.g 2 Sb 2 . The missing entropy above the MFT curve 
in Fig. 10(a) at T = Tn is due to the entropy associated 
with short-range magnetic ordering above Tn- 

We also measured C mag {T) of EuCui. 82 Sb 2 in various 
magnetic fields H applied along the c axis as shown in 



T (K) 


FIG. 11: (Color online) Heat capacity C p of an EuCui. 82 St >2 
single crystal as a function of temperature T measured in 
different magnetic fields H applied along the c axis. 

Fig. 11. Consistent with the x _1 (T) measurements ver¬ 
sus 77 in the upper insets of Fig. 4, the Tn decreases with 
increasing 77. For example, from Fig. 11(e) at 77 = 2.0 T, 
the Tn occurs at 3.8 K which as expected for an AFM 
transition is lower than the zero-field Tn = 5.1 K. It is 
also seen that the heat capacity jump at Tn decreases 
with increase in 77 indicating a second order phase tran¬ 
sition. The 77 — T phase diagram determined from the 77 
dependence of Tn from C P (T) measurements in different 
77 is shown in Fig. 12, where the solid red curve is a fit 

of the function 77 = Hq(1 — ^-) 1//2 to the data. The 
extrapolated critical field at T = 0 is ~ 3.5 T. 


D. Electrical Resistivity 

The in-plane p of EuCui. 82 Sb 2 as a function of T mea¬ 
sured in 77 = 0 are shown in Fig. 13. The magnitude 
and T dependence of p indicate metallic behavior. The 
p decreases almost linearly with decreasing T down to 
20 K below which it tends to be constant, and eventually 
meets an AFM transition at T = 5.1 K as can be seen 
from the inset of Fig. 13, leading to a rapid decrease in 
p below Tn due to loss of spin-disorder scattering. For 
our crystal we find the residual resistivity at T = 1.8 K 
to be po = 23.6 pQcm with a residual resistivity ratio 
RRR = p(300K)//9(1.8K) « 2. 

We analyzed the normal-state p(T) data using the 
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FIG. 12: (Color online) Magnetic phase diagram in the H — T 
plane for EuCui. 82 Sb 2 as determined from the C P (H,T) data 
for a single crystal in Fig. 11, where H is the internal field 
along the c axis. The solid red curve is a fit of the empirical 
function H = -Ho(l — r^-) 1 ^ 2 the data. 



FIG. 13: (Color online) In-plane electrical resistivity p of an 
EuCui. 8 2 Sb 2 single crystal as a function of temperature T 
measured in zero magnetic field. The solid curve represents 
the fit by the Bloch-Griineisen model in Eqs. (38) for 6 K < 
T < 300 K. Inset: Expanded plot of the low-T p data. 


Bloch-Griineisen model for the resistivity arising from 
scattering of electrons from acoustic phonons, 60 

/ T \ 5 / ,0 r/ t t-5 

„ bo (T) = 4K(0 R ) (—) l T — w — )dx , 

(38a) 

where 0r is the Debye temperature determined from re¬ 
sistivity data and 7?.(0r) is a material-dependent con¬ 
stant. The p(T) data for EuCui. 82 Sb 2 were fitted by 

(38b) 


where (po + p s d) is the sum of the residual resistivity po 
due to static defects in the crystal lattice and the spin- 
disorder resistivity p s d due to the presence of disordered 
magnetic moments. The function f(y) of y = T/0r is 
defined by 38 ’ 52 


f(y) = 


Pbg(T) 


Pbg(T = 0r) 


r-l/y 


(38c) 


= 4.226 259 y 5 


■ dx 


I o (e x - 1)(1 - e~x) 

where at T = 0r the electrical resistivity p(0r) is 

Pbg(T = 0 R ) = 0.9 464 635 ft(0 R ). (38d) 


Thus to fit the p(T) data requires three indepen¬ 
dent fitting parameters (po + Psd)) p(©r) and 0r us¬ 
ing Eqs. (38b) and (38c). The solid curve in Fig. 13 
shows the fit of the p(T) data by Eqs. (38b) and (38c) 
for 6 K < T < 300 K using the analytic Pade approx- 
imant fitting function for f(y). 52 The fitted parameters 
are po + p s d = 24.26(2) pQ cm, p(0r) = 10.2(2) pDcm, 
and ©r = 143(2) K. The 7?.(0r) can be calculated from 
the value of p(0r) using Eq. (38d), yielding 1Z(Qr) = 
10.8 pfl cm. Then p s d is calculated from the value of 
Po + Psd using the above 1.8 K value po = 23.6 pfl cm 
which gives p s d ~ 0.7 pfl cm. 


V. PHYSICAL PROPERTIES OF EuCu 2 As 2 
CRYSTALS 

A. Magnetization and Magnetic Susceptibility 

The ZFC and FC %(T) data for EuCu 2 As 2 single crys¬ 
tal measured at different H applied along the c axis 
(Xo H || c) and in the ab plane (xab, H _L c) are shown 
in Figs. 14 and 15. As can be seen from Fig. 14, at low 
applied H = 0.01 T, the %(T) data exhibit well pro¬ 
nounced AFM ordering features at 7 n = 17.5 K for both 
AT || c and H _L c without any hysteresis in the ZFC and 
FC x(T) data. Like EuCui. 82 Sb 2 , the Tn of EuCu 2 As 2 
decreases with increasing H as shown in the insets of 
Fig. 15. Our Tn = 17.5 K is significantly larger than the 
value of 15 K previously reported by Sengupta et al. 40 
At H = 0.1 T we observe a weak change in slope of y(T) 
near 12 K [see inset of Fig. 15(a)] for H _L c, the origin 
of which is not clear. 

The Xab(T) of EuCu 2 As 2 in Fig. 14 exhibits an anoma¬ 
lous strong increase for Tn < T < 20 K that may be due 
to the buildup of FM correlations beyond MFT. This 
behavior is not predicted by MFT and is in contrast to 
the data above Tn for EuCui. 82 Sb 2 in Fig. 3. It is seen in 
Fig. 14 that Xab/Xc « 1.7 at T « Tn. However, at lower 
T < 13 K, Xab < Xc which indicates that the easy axis or 
easy plane is perpendicular to the c axis. This easy-plane 
behavior in EuCu 2 As 2 is the same as observed above in 
EuCui.82Sb 2 . 


p(T) = (p 0 + p sd ) + p(0r)/(T/ 0 r ) 
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FIG. 14: (Color online) Zero-field-cooled (ZFC) and field- 
cooled (FC) magnetic susceptibility x of an EUCU 2 AS 2 single 
crystal as a function of temperature T in the temperature 
range 1.8-30 K measured in a magnetic field H = 0.01 T 
applied along the c axis (y c , H || c) and in the ab plane 

(Xab,H J_ c). 


The x(T) data in the PM state above ~ 50 K are 
well represented by the Curie-Weiss law (1). Linear 
fits of x -1 (T) measured at H = 3 T in the tempera¬ 
ture range 60 K < T < 350 K are shown by straight 
red lines in Fig. 15. The fitted parameters are C = 
7.44(1) cm 3 K/mol and 6 p b = +19.0(2) K for Xab and 
C = 7.64(2) cm 3 K/mol and 9 £ = +17.2(5) K for \c- 
The effective moments obtained from the values of C are 
/i e ff = 7.72(1) /r B from Xab and g eS = 7.82(2) from 
Xc- Again these values of fi e g are similar to the theoret¬ 
ical value of 7.94for a free Eu +2 ion with S = 7/2 
and g = 2, indicating the presence of divalent Eu ions in 
EuCu 2 As 2 . According to Eq. (23), the strongly positive 
values of 9 P ~ Tn indicate the dominance of (negative) 
FM interactions Jij over (positive) AFM interactions be¬ 
tween the Eu spins in EuCu 2 As 2 . 

The isothermal M(H) data of the EuCu 2 As 2 crystal 
measured at different temperatures between 1.8 K and 
300 K for H applied along the c axis (M c , H || c) and in 
the ab plane (M a b, H + c) are shown in Fig. 16. As T 
decreases below 20 K, strong negative curvature occurs 
in M(H) for both H || c and H + c, again indicating the 
presence of dominant FM interactions and correlations. 
It is striking, however, that even at the lowest T = 1.8 K 
this strong negative curvature persists, in contrast to the 
prediction of MFT in Fig. 8 for the perpendicular mag¬ 
netization for an AFM which should be linear in field 
for T < Xn up to the critical field which from Fig. 16 
is of order 2 T. This strong divergence from the predic¬ 
tion of MFT for a planar AFM suggests that the AFM 
magnetic structure is both noncollinear and noncoplanar. 
The noncoplanar component likely consists of canting of 
the Eu moments out of the ab plane in an alternating 



FIG. 15: (Color online) Zero-field-cooled inverse magnetic 
susceptibility x _1 of an EUCU 2 AS 2 single crystal (crystal 1) 
versus temperature T in the T range 1.8-350 K measured 
in a magnetic field H = 3 T applied (a) in the ab plane 
(Xab,H + c) and, (b) along the c axis ( Xc,H || c). The 
insets in (a) and (b) show low-T \ data at different H values. 


way such that the overall magnetic structure is AFM. 

At H = 5.5 T, the saturation moments are = 
6.76 ^ B /f.u. and = 6.95 /t B /f.u., close to the expected 
value of g sa t = gSg B with g = 2 and S = 1/2. The 
M(H) isotherm at 1.8 K is found to exhibit no mag¬ 
netic field hysteresis. The M(H) isotherms at 5 K are 
almost the same as those at 1.8 K. With increasing T, 
/i sa t monotonically decreases as one approaches Tn as ex¬ 
pected. 

The M{H) curve at 1.8 K for H + c is shown sep¬ 
arately in Fig. 17 to illustrate a spin flop transition at 
H ss 0.2 T which is more clearly seen in the derivative 
plots ( dM/dH vs. H) shown in the inset of Fig. 17. The 
slope changes in the M{H) curves for H + c are clearly 
illustrated from the dM/dH vs. H plots at 1.8, 5 and 
10 K in the inset of Fig. 17. In contrast, no such behav¬ 
ior is observed in dM/dH vs. H plot for H || c. This 
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FIG. 16: (Color online) Isothermal magnetization M of an 
EuCu 2 As 2 single crystal as a function of internal magnetic 
field H measured at the indicated temperatures for H (a) 
in the ab plane (M a b,H _L c) and, (b) along the c axis 
(M C ,H || c). 

behavior is consistent with the above y(T) and M(H) 
data that indicate that the ordered Eu +2 moments are 
aligned in the ab plane apart from the above suggested 
canting in and out of that plane. 


B. Heat Capacity 

The C P (T) data for EUCU 2 AS 2 are shown in Fig. 18. 
The C P (T) exhibits a distinct second-order transition 
at 17.5 K as shown in the inset of Fig. 18 due to the 
AFM ordering. The C p attains a saturation value of 
~ 124.5 J/molK at room temperature which is close 
to the expected classical Dulong-Petit value of Cy = 
124.7 J/molK. The C P (T) data were analyzed using 
the Debye lattice heat capacity model in the tempera¬ 
ture range 30 < T < 300 K. The fit of the C P (T ) data 



FIG. 17: (Color online) Isothermal magnetization M of an 
EuCu 2 As 2 single crystal as a function of internal magnetic 
field H measured at 1.8 K for H JL c. Inset: The derivative 
dM/dH vs. H plot at 1.8, 5 and 10 K for Tie and H || c. 



FIG. 18: (Color online) Heat capacity C p of an EuCu 2 As 2 
single crystal as a function of temperature T in the T range 
1.8-300 K measured in zero magnetic field. The solid curve is 
the sum of the contributions from the Debye lattice heat ca¬ 
pacity CV Debye (T) and electronic heat capacity yT according 
to Eq. (35a). Inset: Expanded view of the low-T C P (T) data 
and the estimated lattice contribution of SrCu 2 As 2 (Ref. 38) 
before and after correcting for the difference in formula masses 
of EuCu 2 As 2 and SrCu 2 As 2 . 


by Eqs. (35) using the analytic Pade approximant fitting 
function 52 is shown by the solid curve in Fig. 18 where 
again we used 7 and ©d as adjustable parameters for the 
fit with 30 K < T < 300 K, yielding 7 = 4(3) mJ/mol K 2 
and ©d = 241(3) K. 

The C mag (T) and 5 , mag (T) data for EUCU 2 AS 2 are 
shown in Fig. 19 which were obtained after subtract¬ 
ing the lattice contribution to C P {T ) as done above for 
EuCui. 8 2 Sb 2 . The mass-normalized lattice contribution 
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FIG. 19: (Color online) (a) Magnetic contribution to the heat 
capacity C mag for EUCU 2 AS 2 plotted as C' ma g(T)/T versus 
T. The solid curve represents the MFT prediction of Cmag 
for S = 7/2 and Tn = 17.5 K. (b) Magnetic contribution 
S m a,g(T) to the entropy. 


for EUCU 2 AS 2 is shown in the inset of Fig. 18. Here again 
we observe from Fig. 19(a) that C mag (T) is nonzero above 
the In up to 30 K due to the presence of magnetic correla¬ 
tions. The appearance of the plateau in Cmag (T < Tn)/T 
is consistent with MFT. 22 The C mag (T)/T calculated for 
In = 17.5 K using MFT is shown as the solid red curve in 
Fig. 19(a). Here again the presence of magnetic fluctua¬ 
tions well above the Tn is manifested as a reduced value of 
C , ma g (T)/T below Tn compared to the calculated MFT 
behavior. The S mag (T ) data shown in Fig. 19(b) were 
determined by integrating C mag (T)/T versus T, and is 
found to attain a value of 15.6 J/molK (« 0.901? hi 8) 
at 22 K and 16.1 J/molK (« 0.93I?ln8) at 30 K, thus 
confirming the divalent state of Eu atoms in EUCU 2 AS 2 . 

The C'p(T) data of EUCU 2 AS 2 measured under the ap¬ 
plication of different H along the c axis are shown in 
Fig. 20. It is seen that Tn decreases as H increases. For 
example at H = 1.0 T, the Tn = 13.5 K compared with 



FIG. 20: Heat capacity C p of an EUCU 2 AS 2 single crystal as 
a function of temperature T measured in different magnetic 
fields H applied along the c axis. 



T(K) 


FIG. 21: Magnetic phase diagram in the H — T plane for 
EUCU 2 AS 2 as determined from the C P (H,T) data in Fig. 20, 
where H is the internal field along the c axis. The solid red 
curve is a fit of the empirical function H = Ho(l — 2^-)°' 81 to 
the data. 

Tn = 17.5 K at H = 0. At the same time the heat 
capacity anomaly becomes broader with increasing H. 
We also observe that with increasing H a broad peak 
appears near 20 K for H = 1.0 T [Fig. 20(d)] which 
we suggest is due to field-induced FM correlations above 
Tn(7I). At H = 3.0 T no clear transition is observed 
in the C p (T) data. The H — T magnetic phase diagram 
determined from the H dependence of Tn obtained from 
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FIG. 22: (Color online) In-plane electrical resistivity p of an 
EuCu 2 As 2 single crystal versus temperature T in H = 0. The 
solid curve is the fit by the Bloch-Griineisen model in Eqs. (38) 
for 18 K < T < 300 K. The dashed curve is an extrapolation 
of the fit to T = 0. 


the C p (H, T) measurements in Fig. 20 is shown in Fig. 21. 
The critical field at which AFM is destroyed at T = 0 K 
is estimated from the extrapolated red curve in the fig¬ 
ure to be r; 2.6 T. In contrast, for the polycrystalline 
EUCU 2 AS 2 sample in Ref. 41 with Tn « 15 K the extrap¬ 
olated critical field is ~ 1.8 T. The differences between 
the 7 n and critical field values in Refs. 40 and 41 and 
ours are correlated with the fact that their samples were 
polycrystals whereas ours are single crystals. The differ¬ 
ent preparation conditions evidently lead to differences 
in compositions and/or defect concentrations, resulting 
in differences in the physical properties. 


C. Electrical Resistivity 

The in-plane p(T) data of EuCu 2 As 2 measured in zero 
magnetic field are shown in Fig. 22. The data exhibit 
metallic behavior with po = 14.0(1) gQcm at T = 1.8 K 
and RRR ~ 3.5. As shown in the inset of Fig. 22, a 
sharp transition is seen at Tn = 17.0 K in the p{T) data 
which decreases at lower T due to a reduction in spin- 
disorder scattering. The p(T ) data are well represented 
by the Bloch-Griineisen model. The fit of the p(T) data 
by Eqs. (38) is shown by solid curve in Fig. 22 in the 
temperature range 18 K < T < 300 K using the analytic 
Pade approximant fitting function for pbg(T ). 52 The fit¬ 
ting parameters are po + p s d = 17.44(4) pf 1 cm, p(0r) = 
21.9(3) pflcm, and 0r = 223(3) K. We obtained 7?.(0r) 
= 23.1 pflcm using Eq. (38d). The p s d ~ 3.4 pficm is 
obtained from po+Psd using p o = 14.0(1) pfl cm at 1.8 K. 


VI. SUMMARY 

We have investigated the crystallographic, magnetic, 
thermal and transport properties of EuCu 2 As 2 and 
EuCui. 82 Sb 2 single crystals. The Rietveld refinements 
of x-ray powder diffraction data for crushed crystals in¬ 
dicate full occupancy of the atomic sites in EuCu 2 As 2 , 
whereas Cu vacancies are found in EuCu 2 Sb 2 correspond¬ 
ing to the actual composition EuCui.8 2 Sb 2 . On the other 
hand, no evidence for Cu vacancies in a polycrystalline 
sample of EuCu 2 Sb 2 were found. 49 This difference evi¬ 
dently reflects differences in the sample type (powder ver¬ 
sus single crystal) and preparation conditions. A signifi¬ 
cant concentration of Cu vacancies was previously found 
in the similar compound CaCui. 7 As 2 . 61 

The T dependences of p and C p and the magnitude 
of p indicate metallic ground states in both compounds. 
The p(T > Tn) and C P (T > Tn) data were successfully 
analyzed by the Bloch-Griineisen model and the Debye 
model of lattice heat capacity, respectively. The x(T), 
C P (T), and p(T) demonstrate the occurrence of long- 
range AFM order below Tn = 5.1 K in EuCui. 82 Sb 2 and 
Tn = 17.5 K in EUCU 2 AS 2 . The spin S = 7/2 state of 
Eu atoms in these compounds was indicated both from 
the Curie constant in the Curie-Weiss behaviors of x(^) 
for T > Tn and also by the high-T limit of the mag¬ 
netic entropy S ma . g [T). The x{T) and M(H) data reveal 
anisotropic magnetic properties for both EuCui. 8 2 Sb 2 
and EUC 112 AS 2 at T < Tn- 151 Eu Mossbauer measure¬ 
ments for EuCu 2 Sb 2 at 2.09 K confirmed that the Eu 
spins have a low-T ordered moment of 7 p, b, consistent 
with the value p, sa t = gSpB expected from S = 7/2 and 
9 = 2 49 

The anisotropic magnetic susceptibilities Xab and Xc 
for EuCu 2 Sbi .82 at T < Tn were modeled using molec¬ 
ular field theory. Two AFM structures are equally con¬ 
sistent with the data. The first is a collinear A-type 
AFM structure where the Eu ordered moments in a given 
ab plane are ferromagnetically aligned, and the moments 
in neighboring layers are antiferromagnetically aligned. 
In the second model the Eu ordered moments form a 
noncollinear planar helix with the helix axis being the 
c axis. In both structures the ordered moments are 
aligned within the ab plane. Recent neutron diffraction 
measurements of the magnetic structure of EuCii 2 Sb 2 
showed that the magnetic structure is the collinear A- 
type structure with the ordered moments aligned in the 
ab plane. 49 In Sec. IV B 3, it was shown that the mag¬ 
netic dipole interaction favors the Eu ordered moments 
in the A-type AFM state to lie in the ab plane, consistent 
with the neutron diffraction results and indicating that 
the magnetic dipole interaction is responsible for, or at 
least contributes to, determining the easy axis of the or¬ 
dered moments. On the other hand, whereas EuFe 2 As 2 
shows the same A-type AFM structure 20 of the Eu spins 
as in EuCui. 82 Sb 2 , Co-doped, 62 P-doped, 63 and Ir-doped 
(Refs. 64, 65) EuFe 2 As 2 exhibit FM order of the Eu mo¬ 
ments with the moments oriented along the c axis. From 
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Table IV, FM ordering is not energetically favored com¬ 
pared to AFM ordering by magnetic dipole ordering and 
the ordered moment direction along the c axis is not con¬ 
sistent with the easy axis predicted by this interaction. 
Hence both the FM structure and the easy c axis in these 
doped EuFe 2 As 2 compounds are evidently due to RKKY 
interactions between the Eu local moments. 64 The com¬ 
pound EuZn 2 Sb 2 was reported to exhibit an AFM tran¬ 
sition of the Eu spins at Tn = 13 K. 66 

The magnetic properties of EuCu 2 As 2 at T < Tn 
cannot be understood within molecular field theory for 
collinear and planar noncollinear AFM structures and in¬ 
stead suggest both a noncollinear and noncoplanar AFM 
ground state. This type of AFM state is indicated by the 
strong negative curvature in M(H) isotherms starting at 
small H for both H || c and FI c in Figs. 16 and 17. 

The compound EuPd 2 Sb 2 crystallizes in the 
CaBe 2 Ge 2 -type structure, exhibits an AFM transi¬ 
tion below 6.0 K, a spin reorientation transition near 
4.5 K, and M(H) isotherms with H 1 c at T = 1.8 K, 
but no such field-induced transitions were observed 
for H || c. 43 With what we have learned about the 
predictions of MFT for the anisotropic y(T) of planar 
noncollinear Heisenberg AFMs 22 ~ 24 since Ref. 43 was 


published in 2010, the x(T) data in Fig. 3(a) of Ref. 43 
for EuPd 2 Sb 2 strongly suggest a noncollinear AFM 
ground state such as a planar helical or cycloidal AFM 
structure. Furthermore, from neutron diffraction mea¬ 
surements EuCo 2 P 2 with the ThCr 2 Si 2 -type structure 
exhibits a (planar) helical structure with the Eu ordered 
moments aligned in the ab plane with the helix axis 
being the c axis, where the Co atoms were deduced to be 
nonmagnetic. 6 ' A comparative study of the reason for 
the divergent AFM structures of EuCu 2 As 2 , EuCu 2 Sb 2 , 
EuPd 2 Sb 2 and EuCo 2 P 2 using electronic structure 
calculations 64 would be enlightening. 
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